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We investigate QCD at large p/T by using Za-symmetric SU (3) gauge theory, where p is the quark-number 
chemical potential and T is temperature. We impose the flavor-dependent twist boundary condition on quarks 
in QCD. This QCD-like theory has the twist angle 6^ as a parameter, and agrees with QCD when 9 = 0 and 
becomes symmetric when 0 = 27r/3. For both QCD and the Zs-symmetric SU{3) gauge theory, the phase 
diagram is drawn in p-T plane with the Polyakov-loop extended Nambu-Jona-Lasinio model. In the Z 3 - 
symmetric SU (3) gauge theory, the Polyakov loop p is zero in the confined phase appearing at T < 200 MeV. 

The perfectly confined phase never coexists with the color superconducting (CSC) phase, since finite diquark 
condensate in the CSC phase breaks Z 3 symmetry and then makes p finite. When p > 300 MeV, the CSC phase 
is more stable than the perfectly confined phase at T < 100 MeV. Meanwhile, the chiral symmetry can be broken 
in the perfectly confined phase, since the chiral condensate is Z 3 invariant. Consequently, the perfectly confined 
phase is divided into the perfectly confined phase without chiral symmetry restoration in a region of p< 300 
MeV and T < 200 MeV and the perfectly confined phase with chiral symmetry restoration in a region of p> 300 
MeV and 100 < T < 200 MeV. The basic phase structure of Zs-symmetric QCD-like theory remains in QCD. 

We show that in the perfectly confined phase the sign problem becomes less serious because of = 0, using 
the heavy quark theory. We discuss a lattice QCD framework to evaluate observables at 6^ = 0 from those at 
9 = 27r/3. 

PACS numbers: ll.30.Rd, 12.40.-y 


I. INTRODUCTION 


Quantum chromodynamics (QCD) has a lot of interesting 
phenomena particularly at small temperature (T) and large 
quark-number chemical potential (p). The phenomena may 
affect the structure of neutron stars in its inner core. The 
color superconducting (CSC) phase that appears only at large 
p/T is a typical example. However, lattice QCD (LQCD) 
simulations as the first-principle calculation have the well- 
known sign problem, since the fermion determinant becomes 
complex for real p and this makes the importance-sampling 
method unfeasible in Monte Carlo simulations. Therefore, 
several methods such as the reweighting method 10], the Tay¬ 
lor expansion method00 and the analytic continuation from 
imaginary to real p were proposed so far to circumvent 
this problem, but these methods are considered to be reliable 
only at p/T < 1. Recently, remarkable progress has been made 
with the complex Langevin d yna mics llOt - tl^ and the Lef- 
schetz thimbles approach OfN[l^ . However, the results are 
still far from perfection. Hence the effective model approach 
becomes useful. Actually, QCD at large p/T has been pre¬ 
dicted with effective models such as the Polyakov loop ex¬ 
tended Nambu-Jona-Lasinio (PNJL) model iflA - l^ . 
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The LQCD and effective model approaches have a funda¬ 
mental problem on the definition of quark confinement, as 
mentioned below. It is widely believed that the confinement 
mechanism in SU {N) gauge theory is governed by sym¬ 
metry and confinement and deconfinement phases are defined 
by the Polyakov loop p ll^ as an order parameter for Zn 
symmetry. This is true for pure gauge theory, since Zm sym¬ 
metry is exact. The symmetry is preserved at small T, but 
spontaneously broken at high T. This makes the confinement- 
deconfinement transition first-order in the case of W > 3. In 
SU{N) gauge theory with fundamental fermions, Zm sym¬ 
metry is not exact any more, so that the definition of confine¬ 
ment and deconfinement phases becomes obscure. In fact, in 
QCD with N = 3, p is always finite at T > 0 and the decon¬ 
finement transition is crossover for zero p m. 

An attempt to answer this problem was made recently by 
modifying the fermion boundary condition in Refs. ll32l - [^ . 
Three degenerate flavor QCD was extended by imposing the 
flavor-dependent twist boundary condition (FTBC) 

=/3,x) = —= 0,x) (1) 

on quarks, where /3 = 1/T and Ty is the quark held with 
flavor /. When the twist angles 6 *y are set to 

01 = 0, 02 = -0, 03 = 0, (2) 

the QCD-like theory agrees with QCD when 0 = 0, and be¬ 
comes Z 3 symmetric when 0 = 27r/3, since Of is transformed 
into Of-I by the Z 3 transformation but 6 */-i returns to Of 
by relabeling /. The Z 3 -symmetric QCD-like theory with 
0 = 27r/3 is referred to as Z 3 -QCD in this paper. 
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Confinement and deconfinement phases are clearly defined 
with tf in Z 3 -QCD. The QCD phase diagram can be under¬ 
stood as a remnant of the Z 3 -QCD phase diagram. In Z 3 - 
QCD, it is obvious that there exists a confinement phase at 
small T and a deconfinement phase at high T, because </? = 0 
in the low-T limit and 1 in the high-T limit. We call the phase 
where </? = 0 perfectly confined phase in this paper. The de¬ 
confinement phase transition is either first-order or second- 
order. The order of the phase transition was investigated by 
applying the FTBC to the PNJL model, and was found to be 
first-order at /i = 0 ll32l - [^ . LQCD simulations for 2 h- 1 flavor 
quarks show that is quite small at zero p, and low T 
The confinement property in 2 h- 1 flavor QCD is considered to 
be a remnant of the perfectly confined phase in Z 3 -QCD. 

Z 3 -QCD is useful also to investigate the relation between 
confinement and other mechanism such as chiral symmetry 
breaking and CSC. Particularly for the relation between con¬ 
finement and CSC, we can make clear discussion by using 
symmetry. The diquark condensate A as an order parameter 
for CSC is not Z 3 invariant. This means that Z\ is an order 
parameter for both CSC and Z 3 symmetry. The diquark con¬ 
densate is then zero whenever p is zero, and hence the CSC 
phase with finite A never coexists with the perfectly confined 
phase with zero p. The correlation between p and the diquark 
condensate is thus clearly understood in Z 3 -QCD. The cor¬ 
relation in QCD can be understood as a remnant of that in 
^3 -QCD. Meanwhile, the chiral condensate as an order pa¬ 
rameter for chiral symmetry is invariant. Hence, a chiral 
transition can take place in the perfectly confined phase. In 
fact, the PNJL model shows that this really happens at finite p 
ll3^ . In the PNJL analysis, however, diquark effects were not 
considered. 

In this paper, we try to understand QCD at high p/T from 
Z 3 -QCD. The phase diagram is drawn in p-T plane for both 
QCD and Z 3 -QCD by using the PNJL model with the FTBC. 
In the case of Z 3 -QCD, a perfectly confined phase with p = Q 
and A = 0 appears at T < 200 MeV, and the chiral restora¬ 
tion occurs at p 300 MeV in the phase, as expected. 
Hence, the perfectly confined phase is divided into the per¬ 
fectly confined phase without chiral symmetry restoration in 
a region of p < 300 MeV and T < 200 MeV and the per¬ 
fectly confined phase with chiral symmetry restoration in a 
region of p > 300 MeV and 100 < T < 200 MeV. In a region 
of p > 300 MeV and T < 100 MeV, CSC phases with finite 
A and small p come out, since they cannot coexist with the 
perfectly confined phase and more stable than the perfectly 
confined phase with chiral symmetry restoration. This basic 
structure remains in QCD. Second we show that in the per¬ 
fectly confined phase the sign problem may become less se¬ 
rious because of ip = 0, using the heavy quark model. This 
may make LQCD simulations feasible in the perfectly con¬ 
fined phase. Particularly when the system is in the perfectly 
confined phase without chiral symmetry restoration, the sys¬ 
tem changes continuously from Z 3 -QCD to QCD by varying 
9 from 27r/3 to 0. Using this property, we propose a way of 
estimating observables for QCD from those for Z 3 -QCD. 

This paper is organized as follows. We recapitulate the 
FTBC and Z 3 -QCD in Sec. Ulland the PNJL model in Sec. 


Hm Numerical results are shown by using the PNJL model and 
LQCD simulations in Sec. |IV] The sign problem in Z 3 -QCD 
is discussed in Sec. [V] Section IVlIis devoted to a summary. 

II. Z3-QCD 

In this section, we recapitulate the FTBC and Z 3 -QCD, fol¬ 
lowing Refs, We consider SU{N) gauge theory with 

N degenerate flavor quarks. The Lagrangian density C in Eu¬ 
clidean spacetime is 

^ 1 

^ = E (3) 

/=! ^ 

where m is the current quark mass for all flavors and D^, = 
di, — iAy for the gauge field and is the field strength 
tensor. The temporal boundary conditions for quarks are 

^/(■r =/3,x) == 0,x). (4) 

The Lagrangian density Q is invariant under the Zm (large 
gauge) transformation, 

Ff ^ Ff = UFf, 

^ A'^ = UA^U-^+i{dM)U-\ (5) 

but the boundary condition (IDl is changed into 

Ff{T = P,x.) = -F'^'^^/^Ff{T = 0,x). ( 6 ) 

Here 

U{x,t) = exp(faafa) (7) 

is an element of SU{N) group characterized by real functions 
aa{x, t) satisfying the boundary condition 

U{x,P) = exjp(—i2TTk/N)U{x,0) (8) 

for any integer k. 

Now, we consider the following FTBC instead of (|4|i: 

Ff{T = I3,:x.) =-e~"^fFf{T = 0,x.), (9) 

where 

27T 

0/=0i + ^(/-l) (/ = 1,2,...,A^/), (10) 

with Nf = N. Under the Zpf transformation 0 , the bound¬ 
ary condition 0 is changed into 

!f/(T =/3,x) == 0 ,x), ( 11 ) 

where 

Stt 

e'f=9, + —if-l-k) (/ = l,2,...,iV). (12) 

The boundary condition (fTST i returns to the original one (fTOl i 
by relabeling the flavor index f — kasf. Hence SU (TV) gauge 
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theory with the FTBC @ has symmetry exactly. SU(N) 
gauge theory with I x N flavor fundamental quarks also has 
Zm symmetry for any positive integer I, when the FTBC dUl 
with Nf = IN is imposed on the fermions ll3^ . 

When the fermion fields iPf are transformed as ifssi] 

^f ^ exp {—i9fTT)'Ff, (13) 

the boundary condition (|9|l returns to the ordinary one dH, but 
the Lagrangian density C is changed into 

= E (14) 

/=! ^ 

where D® = 9^ — i{Ai, + OSv^rT) and iOT is the flavor- 
dependent imaginary chemical potential given by the matrix 

e = diag(6»i,02, • • • ,9n) 

= diag(6»i,0i-I-27r/iV, • • • , 

01 + 2(/ - • • • , 01 + 27r(iV - l)/N) (15) 

in flavor space. In the case of TV = 3, is nothing but the 
Lagrangian density of Z 3 -QCD. The flavor-dependent imag¬ 
inary chemical potential breaks SU(N) flavor symmetry and 
associated SU{N) chiral symmetry HH]. In the chiral 
limit, global SUy{3) x S'?7a( 3) symmetry is broken down to 
(17(l)v)"®(C/(l)A)"i3l. The symmetry is even broken into 
((7(l)v)^, when chiral symmetry is spontaneously violated. 


III. PNJL MODEL WITH FLAVOR-DEPENDENT TWIST 
BOUNDARY CONDITION 


In this section, we explain the Polyakov-loop extended 
Nambu-Jona-Lasinio (PNJL) model with the FTBC (l9]l and 
diquark effects. Taking the Polyakov-gauge and treating A 4 as 
a background gauge field, one can get the model Lagrangian 
density for three degenerate flavors as 

4njl = E 4 ( 7.4 + H'^/+ 4 }l+Z^>( 16 ) 

f—u,d,s 


where Df, = — iS^^ 4 {A 4 + 9T). In (fThl l. U isa function of 

ip and its complex conjugate p* defined by 


o 

with the condition (/>i -f ^!)2 + <t>p = 0 . 
We take the following form as U: 


(17) 




aiT) , 


-I- 6 (T) ln(l — 6 pp* -f -I- p*^) — S^pp*)"^) 
'Tn\ /Ton2 


(18) 
rn\3 


i(T) = ao + ai(Y) +“2 (y) , &(7’) = &3(y) ■ 


( 19 ) 


The parameters in hi are determined from LQCD data 
in the pure gauge (heavy quark) limit. The Polyakov-loop 
potential has a first-order deconfinement phase transition at 
T = To in the limit, and hence Tq = 270 MeV. In the case 
of finite quark mass, however, the PNJL model with this value 
overestimates the pseudocritical temperature Tc ~ 160 MeV 
at zero p determined by full LQCD lEJSl. We have then 
rescaled Tq to 195 MeV to reproduce Tc ^ 160 MeV j^ . 
The parameters in the Polyakov-loop potential are summa¬ 
rized in TableHIa). 


(a) 

do 

ffli 

02 

63 

ro(MeV) 

3.51 

-2.47 

15.2 

-1.75 

195 

(b) 

m/(MeV) 

A(MeV) 


GdA^ 

KA^ 

5.5 

602.3 

1.835 

!GsA 2 

12.36 


TABLE I: Summary of the parameter set in the PNJL model for 
the case of V = 3. Panels (a) and (b) show the parameters in the 
Polyakov-loop potential and the NIL sector, respectively. 


In m, 4k stands for the effective quark-antiquark and 
quark-quark interactions H: 


/’int 

^NJL 


= -Gs 


a—0 


oc—u,d,s c—r,g,b 

x[(Tc)k75e""'^e„.,(T)“] 


+K 


detT/(l -F 75 )T// -F h.c. 


( 20 ) 


where the Aa are the Gell-Mann matrices in flavor space, Gs 
and G^) are coupling constants of the four-quark interactions 
and K is a coupling constant of the Kobayashi-Maskawa-’t 
Hooft determinant interaction |4^ . The quark field 
carries color (a = r, g, b) and flavor {a = u, d, s) indices. The 
values of coupling constants, current quark mass and three 
dimensional momentum cutoff A are tabulated in Table Ub). 
The coupling constants and the cutoff were determined to re¬ 
produce empirical values of g'- and 7r-meson masses and tt- 
meson decay constant at vacuum when m-u = md = 5.5MeV 
and rris = 140.7MeV In this paper, however, we con¬ 
sider a symmetric current quark mass ofm = TO/ = 5.5 MeV 
to preserve SU(3) flavor symmetry. 

Taking the mean field approximation, one can get the ther¬ 
modynamic potential per volume as 


17 = Qq + U + U 


2NNf 

-E 


d?p 


Tj-F (1-F 


(27r)3 I 

where the mean-field potential part U is given by 

U = Gu Y. |4|4 2 Gs E 

/—1,2,3 f—u,d,s 

-AKauUdCTs 


, ( 21 ) 




( 22 ) 
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with chiral condensates cr/ = (S'/'?/) for / = u,d,s and 
diquark condensates Ai = where no 

sum is taken over I on the right hand side. In (l221 i. the Ej are 
the quark spectra which depend on the absolute value of three- 
dimensional quark momentum p, the effective quark mass 

Mf = ruf — 4:Gscrf -|- 2K(j fcr fn, 

if f", f + f") (23) 

the diquark condensate (multiplied by —200 ) 

Ai = -2GdA, (24) 

and the effective chemical potential 

fi'f = fi + i(j)cT + iOfT. (25) 

It is difficult to obtain the explicit forms of Ej analytically, 
but we can determine the Ej by solving the eigenvalue equa¬ 
tion for Dirac operator numerically. Using these spectra, we 
can calculate il and then find the solutions <Jf,Ai and (p that 
minimize 17. Obviously, cr/ is invariant under the Z3 transfor¬ 
mation ©, but Ai is not. 

IV. NUMERICAL RESULTS 

In this section, we show numerical results of the PNJL 
model. For later convenience, we put 

9^ = di = 9, 9d = 92 = -9, 0 , = 03 = 0 . (26) 

The fermion boundary condition (|2^ agrees with the FTBC 
([Tol l when 9 = 27r/3 and the standard antiperiodic boundary 
condition when 0 = 0 . 

In Z 3 -QCD with 0 = 27r/3, the flavor symmetry is vio¬ 
lated with the FTBC (as the mentioned above about the flavor- 
dependent imaginary chemical potential), but u and d quarks 
are symmetric under the interchange between u and d. Since 
the boundary condition is changed by the Z 3 gauge transfor¬ 
mation, flavor and diquark indices are renamed so that the 
conditions Mi < M 2 < M 3 and |Z\i| < |Z\ 2 | < IZial can 
be satisfied. 

Figure [T] shows T dependence of the absolute value of di¬ 
quark condensates Z\i, Z \2 and As at p, — 340 MeV for 
Z 3 -QCD with 0 = 2tt/3. There appear a variety of CSC 
phases Below T — Ti =10 MeV, all the diquark con¬ 
densates are finite, indicating that it is the color flavor locking 
(CFL) phase. Two of three condensates are finite in the region 
Ti < T < T 2 = 80 MeV, and one of three is finite in the 
region T 3 = 223 < T < T 4 = 419 MeV. This means that the 
former is the uSC phase and the latter is the 2SC phase. 

In Fig. |2] the average values of dynamical quark masses 
and the absolute value of diquark condensates, namely, M — 
(Mi-|-M 2 -|-M 3 )/ 3 andZ\ = (|Z\i|-|-|Zi 2 |-|-| 4 \ 3 |)/ 3 , are plot¬ 
ted as a function of T, together with <P = \(p\, at p = 340 MeV 
for Z 3 -QCD with 0 = 27r/3. There are three equivalent solu¬ 
tions to ip: ip = with (f> = 0, ±27r/3 ( or (/> = tt, ± 7 r/ 3 ). 
The solutions of (j) = ±27r/3 (± 7 r/ 3 ) are Zs images of the 


0.2 
0.18 
0.16 
0.14 
5: 0.12 

S 0.1 

^ 0.08 

0.06 
0.04 
0.02 
0 

0 0.1 0.2 0.3 0.4 0.5 

T (GeV) 

Fig. 1: T dependence of |zii| at /i = 340 MeV in the PNJL model 
with 6 = 27r/3. The diquark condensates Ai, A 2 and A 3 are 
denoted by dashed, solid and long dashed lines, respectively. Be¬ 
low r = 10 MeV, three diquark condensates exist. In a region of 
T = 10 ~ 80 MeV, Z\i = 0 and Z \2 = Zla. 


solution of (^ = 0 (tt). In fact, they are thermodynamically 
equivalent and are transformed from one to another by the Z 3 
transformation. We can then take the solution of (^ = 0 (tt), 
and if is real. As mentioned in Sec. |I] finite diquark con¬ 
densates in CSC phases break Zs symmetry and hence induce 
finite <P. In the CFL and uSC phases below T 2 , ^ is tiny but 
finite, as expected. Thus, CSC phases do not coexist with the 
perfectly confined phase, but with the almost-confined phase. 
Also in the 2SC phase at T 3 < T < T 4 , <P is finite as ex¬ 
pected, but it is large. Thus quarks are deconfined in the 2SC 
phase. In the region T 2 < T < Ts, ‘P is zero and hence di¬ 
quark condensates are also zero. This is the perfectly confined 
phase. Above T 4 , (p is finite but A is zero. This is the pure 
deconfinement phase without diquark condensate. 

Figure [3 is the same as Fig. [T] but for QCD with 0 = 0. 
Below T = Ti =32 MeV, all Ai are finite, indicating that it 
is the CFL phase. In the region Ti < T < T2 = 86 MeV, 
one of three appears. It is the 2SC phase. The phase structure 
is thus much simpler at 0 = 0 than at 0 = 27r/3. Figure|4]is 
the same as Fig. |2l but for 0 = 0. In this case, ip is always 
real and p = \p\ = (p. Since Zs symmetry is not exact in 
this case, (p is always finite at T > 0. Thus we cannot define 
the confinement and deconfinement phases clearly. Since <P is 
small in the CFL and 2SC phases, one can consider that the 
CSC phases are in the almost-confined phase. 

Figure|3shows the phase diagram in p-T plane for the case 
of 0 = 27r/3. A variety of phases are present in this case. 
When T increases from zero with p fixed at a large value, say 
350 MeV, the CFL, uSC, perfectly-confined, 2SC and pure 
deconfinement phases appear. We also see that the chiral tran¬ 
sition takes place in the perfectly confined phase. 

Figure | 6 ] is the same as Figure |3 but for 0 = 0. The 
uSC phase disappears and the 2SC goes down to lower T, 
while at small T the CFL phase remains and the perfectly- 
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Fig. 2: T dependence of M, A, and $ at = 340 MeV in the 
PNJL model with 9 = 2tv/3. Here, M, A and ^ are denoted by 
dashed, solid and long dashed lines, respectively, and M and A are 
normalized by the constituent quark mass Mo at the vacuum. 


Fig. 4: T dependence of M, A and at /r = 340 MeV in the PNJL 
model with 6 = 0. Here M and A are divided by Mq. See Fig. [^for 
the definition of lines. 



Fig. 3: T dependence of |Zii| at fi = 340 MeV in the PNJL model 
with 9 = 0. See Fig. [T]for the definition of lines. 


confined phase becomes an almost-confined phase without di¬ 
quark condensate. The almost-confined phase and the CFL 
phase at low T can thus be regarded as remnants of the per¬ 
fectly confined phase and the CFL phase in Z^-QCD. 


Fig. 5: Phase diagram in the PNJL model with 9 = 27r/3. The 
CFL phase exists below the dotted line. The solid line stands for the 
first-order chiral phase transition; M is large on the left side of the 
solid line, but small on the right side. The dashed double-dotted line 
denotes the first-order deconfinement phase transition; the system is 
in the confinement phase below the line and in the deconfinement 
phase above the line. The perfectly confined phase is labeled by 
“<P = 0, A = 0”, while the CFL, uSC and 2SC phases are by 
“CFL”,“uSC” and “2SC”, respectively. 


following relation : 


V. SIGN PROBLEM IN Z3 -QCD 

It is well known that LQCD has the sign problem at real 
p and it makes the importance sampling method unfeasible. 
Also in the case of Z 3 -QCD, the sign problem exists in prin¬ 
ciple. In fact, the determinant of Dirac operator A4 has the 


[det Ad(/r/)]* = detAd(— 

= det[U — — i0/T)74-F TO/] 

f—u,d,s 

= det[D — + mf] 

f—u,d,s 

= det 


( 27 ) 
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quark mass m. In general, detAi{fJ,f ) is not real, since the ex¬ 
ponential factor with Tr[{7x] is not equal to that with its conju¬ 
gate Tr[{7x]^- For the configurations that satisfy the condition 
Tr[[/x] = 0 , however, only the first and fourth terms remain 
on the right hand side of of (l30l l. and hence det A4(p/) be¬ 
comes real. In Z 3 -QCD, we use fj,f instead of p,, but the factor 
g± 3 /i//T _ g± 3 /i gjjjj when 9 = 27r/3. This indicates 
the possibility that the sign problem may not be serious in the 
perfectly confined phase, since the gauge configurations are 
concentrated on the vicinity of = 0. One may then perform 
the importance sampling in lattice Z 3 -QCD. 

Once a thermodynamical quantity 0{9) and it derivatives 
with respect to 9 are obtained at 6 * = 27r/3, the quantities at 9 
less than 27r/3 are obtainable by using the Taylor expansion 


(GeV) 


Fig. 6: Phase diagram in the PNJL model with 6 — 0. The solid line 
(thin dotted line) stands for the first-order (crossover) chiral transi¬ 
tion. The thin dashed douhle-dotted line denotes the line of$ = 0.5. 
The almost-confined phase is labeled by < 0.5”, while the de- 
confined phase is by > 0.5” and the CFL and 2SC phases are by 
“CFL” and “2SC”, respectively. 


for = ^ + i9fT satisfying the condition (l26T l. where the 
third equality has been obtained by relabeling the /. Hence, 
the determinant detAd(/i/) is not real. The present system 
thus has the sign problem, although the partition function is 
real, since 


ZiiifY = = Z{f,f), (28) 

where the first equality comes from (|27] | and the second one 
from charge conjugation. Note that Eq. (l28T l is true for any 
value of 9. 

Although Z 3 -QCD has the sign problem at real p, there is a 
possibility that it is not serious in the perfectly-confined phase. 
To see this clearly, we consider the heavy quark model Gl 
In the model, the fermion determinant of Dirac operator 
is given in terms of the Polyakov line operator Tr[t7x]: 

det Ad(/r/) = det[l -F he'^^'^Ux]det[l + he~^^'^UY\ (29) 

with 


det[l -F /le'^/'^C/x] = 1 + /ie'^/'^Tr[t/x] 

_F;j2g2M/TTr[t/t] + /i3g3M/T^ 

det[l -F = 1 + /ie-^/^Tr[t/J:] 

(30) 

where (7x = (7(x, 0) for the lattice link variable (7(x, 0) de¬ 
pending on the space coordinate x only. The parameter h 
is defined as h = C^* with the lattice-site number Nt in 
the temporal direction and the hopping parameter C, for Wil¬ 
son fermion and l/( 2 m) for staggered fermion with a current 


0 ( 0 ) = y; 

n—0 


1 5”0(6») 
n! 90" 


e=2Tr/3 


(- 1 ) 


(31) 


First we consider the perfectly confined phase with chiral 
symmetry breaking in Z3-QCD. The phase is located in a re¬ 
gion of T < 200 MeV and p < 300 MeV. In QCD with 0 = 0, 
it is very likely that the region includes the nuclear-matter re¬ 
gion at /r ~ 300 MeV and T = 0. 

Figure [T] shows 0 dependence of tP and the quark number 
density Uq at the point (/i,T) = (260 MeV, 100 MeV). At 
0 = 27r/3, this point is in the perfectly confined phase with 
chiral symmetry breaking. We see that the 0-dependence is 
smooth in a range of 0 = 0 ~ 27r/3. Hence, the approach 
mentioned above may be valid for this point because there 
is no singularity. In particular, increases monotonically as 
0 decreases from 0 = 27r/3. This may suggest that only the 
lower derivative terms are needed in the Taylor expansion dSTT l 
to evaluate <P at smaller 0. On the contrary, Ug has a minimum 
at 0 ~ 1.3. In the Stefan-Boltzmann limit, Ug is given by 

Uq = [fi + i9T)T'^ + in + i9T)^ 

-F(/r - i9T)T'^ + (m - iOT)^ 

= 2 - 3tJ.9^T^) , (32) 

and it decreases monotonically as 0 increases from 0 = 0. As 
easily seen in Eq. (l30l l. however, the 0 term is suppressed by 
strong confinement near 0 = 27r/3. As a consequence of the 
suppression, Ug is considered to have a minimum at 0 ~ 1.3. 
The fact that 0 dependence is not monotonic for Ug means that 
the higher derivative terms are necessary in (|3TI ) to evaluate Ug 
at smaller 0 . 

Figure 0 is the same as Fig. |7] but at the point (/i,T) = 
(300 MeV, 100 MeV). At 0 = 27r/3, this point is in the 
perfectly confined phase with chiral symmetry restoration. We 
see that both and Ug have discontinuities as a function of 0 . 
Hence, the Taylor-expansion approach of (ISTT l does not work 
at this point. Thus, the Taylor-expansion approach of (ISTT i 
may work, if the system is in the perfectly confined phase 
with chiral symmetry breaking in the case of Z^-QCD. 

We remark that, in the numerical calculations for Figs. |7] 
and | 8 ] we assumed that the diquark condensate does not ap¬ 
pear in the whole region of 0 = 0 ~ 27r/3. Since the analytic 
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Fig. 7: 9 dependence of ^ and the quark number density Uq at 
(/r, T) = (260 MeV, 100 MeV) in the PNJL model. The solid 
line stands for while the dashed line corresponds to Uq. The Uq is 
divided by 3no, where no is the normal nuclear matter density. 


Fig. 8: 9 dependence of $ and the quark number density Uq at 
(/r, T) = (300 MeV, 100 MeV) in the PNJL model. See Fig. |7] 
for the definition of lines. 


form of the quark spectrum Ej is not known, we must use nu¬ 
merical differentiation to determine the quark number density 
Uq, when Ai ^ 0. We found that it is very difficult to perform 
such differentiation with high accuracy in the low temperature 
region where T < lOOMeV. Then, we adopt the assumption 
mentioned above and performed the approximated calcula¬ 
tions. This assumption is very natural for p = 260MeV where 
the chiral symmetry is broken due to large quark mass M but 
may not be valid for ^ = SOOMeV where M is small. If the 
diquark condensate appears at some value of 0 in the case of 
(/i, T) = (SOOMeV, lOOMeV), Uq and $ is not analytic there. 
Hence, the conclusion that the Taylor-expansion approach of 
(ISTT i does not work at (/r, T) = (SOOMeV, lOOMeV) is not 
changed. 


VI. SUMMARY 

In summary, we investigate QCD at large ^/T by using the 
SU (S) gauge theory with the flavor-dependent twist boundary 
condition (FTBC). The theory agrees with QCD at 6 * = 0 and 
becomes symmetric at 6 * = 27r/S. In Z 3 -QCD, one can 
make clear discussion on the relation between Za-symmetry 
and diquark and chiral condensates. 

In Z 3 -QCD, there can exist a perfectly confined phase 
where = 0. Since the diquark condensate A is an order 
parameter for both Z 3 symmetry and CSC, A is finite (zero) 
when (p is finite (zero). The perfectly confined phase with 
tP = 0 thus never coexists with the CSC phase with A ^ Q. 
Meanwhile, the chiral condensate as an order parameter for 
the chiral transition is Z 3 invariant. Hence the chiral transi¬ 
tion can take place in the perfectly confined phase. 

The phase diagram was numerically investigated for Z 3 - 
QCD with the PNJL model. The perfectly confined phase 
with ^ = 0 and Z\ = 0 appears at T < 200 MeV, and 


the chiral restoration occurs at p Ri 300 MeV in the phase. 
Hence, the perfectly confined phase is divided into the per¬ 
fectly confined phase without chiral symmetry restoration in 
a region of /r < 300 MeV and T < 200 MeV and the perfectly 
confined phase with chiral symmetry restoration in a region 
of p > 300 MeV and 100 < T < 200 MeV. In a region of 

> 300 MeV and T < 100 MeV, CSC phases with flnite Z\ 
and small P appear, since they are more stable than the per¬ 
fectly confined phase with chiral symmetry restoration. This 
basic structure remains in QCD, although the phase structure 
is much more complicated in Z 3 -QCD than that in QCD. The 
phase diagram for QCD is thus a remnant of that for Z 3 -QCD. 

We also showed that in the perfectly confined phase the sign 
problem may become less serious because of ^ = 0 , using the 
heavy quark model. This may make LQCD simulations feasi¬ 
ble in the perfectly confined phase. Particularly when the sys¬ 
tem is in the perfectly confined phase without chiral symmetry 
restoration, the system changes continuously from Z 3 -QCD 
to QCD by varying 0 from 27r/3 to 0. Using this property, 
we proposed the Taylor-expansion method of OTI) to estimate 
observables for QCD from those for Z 3 -QCD. In the perfectly 
confined phase with chiral symmetry restoration, however, the 
system does not change smoothly as 0 varies from 27r/3 to 0. 
Therefore, the lattice QCD framework proposed above may 
not work. Further study is necessary for this region. 

In this paper, we considered three degenerate flavor QCD 
as the first step. Even in the 2 h- 1 flavor case, there exists a 
point of ^ = 0 in /r-T plane for some 6 near 27r/3 JUt]. 
For example, in the case of infinite nis (the two-flavor case), 
such a point appears at 6 * 7 r/ 2 . Hence, we can expect that 

the present Taylor expansion method is applicable for realistic 
2 h- 1 flavor QCD. 
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